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ABSTRACT 


four eneweomaedualttative manner the kinetics of an 
excimer laser, a simple computer model for calculating the 
species populations in a KrF laser cavity is developed; sub- 
sequently a computer program originally developed at the 
Naval Research Laboratory (NRL) is modified to calculate the 
population of the different electronic configurations of 
excited neon present in a XeF laser. When modified, the NR 
Meestam accounts for 135 reactions and requires 9.5 minutes 
of CPU time using the IBM 360/67. The populations obtained 
are applied to the calculation of the index of refraction 
in the laser cavity. The phase shift is determined per unit 
length for the two laser wavelengths; one laser wavelength is 
non-resonance with neon and the other is at resonance. Neon 
is the dominant specie relative to concentration and within 
the population distribution the neon ground state dominates 
Meee tactor of a million. The calculations show that the 
ground state neon determines the index of refraction; an 
exception occurs if the laser wavelength is near resonance 
to any of the transitions in the 3s + 4p array. As long as 
the laser is operated away from the 3s + 4p resonant wave- 
lengths, the phase shift will be negligibly small resulting 
in satisfactory beam quality; the preceding statement is 


valid only for the influence of neon. 
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The rare-gas monohalides, such as XeF, are simple diatomic 
molecules whose properties and emission spectra were essen- 
tially unknown as recently as five years ago. Now they are 
the active media for gas lasers that could provide overall 
electrical efficiency as high as 10 percent in various wave- 


lengths, as depicted in Table lI. 


Table I. Rare-gas Halide Laser Emission Wavelengths 


Peak wavelengths of most intense band (nm) 


Ne Ar Kos Xe 
F Oe 193 248 50, 353 
Cl b 170 222 308 
Br b ieee 2062 282 
i b b Leger 252 


Shins memeogemonstrated laser oseamllation. 


b 


Predissociates, hence emission is weak or unobservable. 


The overall efficiency as explained by J. J. Ewing and cC. A. 
Beau {l] is 


..a product of a quantum efficiency, an upper 
laser level production efficiency, an extraction 
Sentclency, am efficiency for producing the initial 
excited states by the pump, and an energy coupling 
efficiency which describes the efficiency with which 
energy gets from a wall socket into the gas medium. 
The extraction efficiency is comprised of both 
Spectroscopic andkinetic extraction efficiencies 





depending on the ratio of net gain to loss and the 

ratio of rates of stimulated emission of upper laser 

levels to the quenching and spontaneous loss of 

upper levels. Certain rare gas halide lasers have 

potential for high power and high efficiency because 

they have high efficiency in each of the above 
mentioned elements. Depending on the application 
requirements, one may accept a lower efficiency of 

one of these elements if the compromise allows for 

some simplification of the laser. Typically this 

trade off usually involves coupling and excited state 

production efficiencies. 

Rare-gas monohalides belong to a larger class of mole- 
cules broadly known as excimers. An excimer is an atomic 
or molecular aggregate, composed of two atoms or molecules, 
Mate tS DOoUNnd in itS excited state but is unstable or 
slightly bound in the ground state. 

The term "excimer" is standard chemical nomenclature 
for dimers which are bound in the excited state and free in 
the lower state. In standard nomenclature dimer means two 
of the same thing. Thus, xen™, AT5*, etc. are excimers. 
Species that are heteronuclear such as LiXe* or KrF* also 
can have bound excited states and dissociative lower states. 
The accepted chemical nomenclature for such species is 
"exciplex" short for excited complex. Unfortunately, the 
laser community has not uniformly utilized good scientific 
English and typically calls all species with dissociative 
lower states, and even some like XeF with slightly bound 
lower states, ''excimers". 

The excited state of the excimer species can radiate 


in a broad band that is red shifted from the wavelength of 


the parent atomic excitation. The size of the bandshift 


10 





depends on the depth of the excited state well (typically 
1-5eV) and on the repulsion of the lower state at the energy 
minimum of the excimer. In Figure 1 are schematic potential 


energy curves for an excimer. 
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Figure 1. Potential Energy of Excimers (Reproduced from 
Ewame [247 


Excimer emission spectra are useful for lasers because 
the lower level rapidly dissociates. Since the dissociation 
time, about li 4 second, is much less than upper level 


7) 6) TO 8 second, population inver- 


radiative lifetimes, 10 
sions can be readily produced and maintained. 

The excimer laser or exciplex laser offers features 
which are essential for military applications [2-4]. These 
features are as follows: first, the potential for high 


efficiency; second, scalability to large volume and high 


power; and, third, transmission through the atmosphere 


ieee 





using either directly-generated photons or frequency-con- 
verted photons [5]. As a result of the attractive character- 
istics, considerable effort has been and is being directed 
toward identification of key technical issues and establish- 
ment of feasible solutions to the problems identified [6-11]. 
In addition to research directed toward the laser device 
itself, work is also being funded on related and relevant 
aspects, including e-beam technology, UV optics, and time 
dependent laser aerodynamics. For example, see Bennett, 
Seeel. [12] for work on UV optics. 

Systems studies have been conducted by Rocketdyne. Among 
the various aspects of large excimer lasers considered in the 
Rocketdyne reports was beam quality as influenced by non- 
UmELOrmity of index of refraction. The adverse impact of 
poor beam quality upon the overall system was stressed. 

For military applications the optimum pulse duration 
should be a few microseconds [13-15]. The pulse duration 
for some military applications is in contrast to the desired 
pulse duration of a few nacoseconds for inertial confinement 
to achieve thermonuclear fusion [1]. Even if a pulse dura- 
meon Of 10 microseconds is used, the motion of the waves, 
e.g., cathode shock wave and anode wave, in the gas will be 


very small, i.e., of the order of 
2 = at = (300 m/sec) (10° sec) = 3mm (1) 


where a is the acoustic wave speed and t is the vulse dura- 


tion. However, this fact does not suggest that a pulsed 


eZ 





excimer laser will be free of inhomogeneity problems; exten- 
sive experience with pulsed electrical lasers operating in 
the IR indicates that beam quality may be a problem. To 
achieve high average power, a large pulse repetition fre- 
quency, PRF, must be used which requires a flowing laser. 
The higher the PRF, the faster the gas flow is. 

Gas flow in a high power, high PRF, excimer laser will 
probably be subsonic. In subsonic flow [16-18], disturbances 
can propagate upstream. Even if the upstream base flow is 
perfect, the gas moving out of the cavity after a pulse can 
cause An/n within the cavity; n is index of refraction, and 
An is spatial variation of index of refraction. A large 
flush factor eases the problem of a nonuniform base flow; 
however, high efficiency demands a low flush factor. 

Nonuniformity of the lasing medium may have adverse 
effects distinct from the deterioration of beam quality. 

The work of Nerheim and Chen [19] at JPL indicates alts 5 S2cisotte 
levels of performance depending on homogeneity of the gas. 
Pulses into a homogeneous medium resulted in one level of 
power output; pulses into an inhomogeneous medium yielded 
another level (lower naturally) of performance. Details of 
the physics involved in this observation were not unraveled. 
The work by Nerheim and Chen was with an atomic copper laser. 
Consequently the results are not directly transferable to 
rare gas-halide lasers. 

Despite the many unsolved problems regarding excimer 


lasers, its use as a potential weapon constitutes a 


3 





revolution greater than the leap from naval guns to guided 
missiles. Some of the conceivable missions for the excimer 
laser are ballistic missile defense (both ground-based silo 
defense and space-based destruction of ICBMs in the boost 
phase), space-based antisatellite defense, and antiship 
missile defense (ASMD) [20]. The latter has the most pro- 
mise of relatively near-term development. Also, the radia- 
tion from certain excimer lasers has good propagation in 
Pea water. 

Briefly, the laser offers the advantage of depositing 
an extremely high energy on a small area of a target. This 
destructive energy is delivered at the speed of light (399 
million meters per second which is ideal for a defensive 
system against missiles); whereas, a conventional missile 
travels to its target at a speed of a few times the speed 
of sound (roughly 1500 meters per second). 

As alluded to earlier, beam quality has an important 
impact upon the overall laser system. Improvements in beam 
quality can be used to decrease overall laser size, shorten 
time for melting a metal, or decrease laser run time to 
perform a certain task. A well know fact is that beam 
quality depends on favorable illumination of the exit 
aperture; favorable illumination implies near constant 
phase and, to a lesser degree, near constant amplitude. 
Variation of the index of refraction within the cavity 


causes spatial variation in phase at the laser exit. This 
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nonuniform phase causes the peak axis intensity, Tee, Ee 


decrease according to 
E 2 
a = exp [-(27d 5/4) 4) (2) 


where A is wavelength. The rms phase, ¢ is obtained by 


rms ’ 
squaring and averaging the phase shifts across the aperture. 
Equation (2) is the Strehle Ratio [21]. The spatial phase 


fmeratilon 1s 


6 = f nds (3) 


where ds is an incremental distance along a ray in the 
laser cavity. The initial efforts to limit phase distortion 
considered the index of refraction in terms of the Gladstone- 


fee LOrmulation 


Sa ss ees (4) 


where 9 is gas density, and G is the Gladstone-Dale constant; 
hence there is concern for flow and acoustic effects. Hogge 
and Crow [22] state that flow uniformity and reduced acoustic 


effects have achieved Ap/o = on 


, which is the required 
meow Uniformity. 

It has been suggested, however, that the requirement for 
flow uniformity should be stated in terms of An/n due to 


Mem@uoLrOrm pumping in the cavity [23]. In a typical flow 


Mixture, the constituents of the medium (Ne :Xe: NF.) are 


LD 





present in a fixed ratio (750:3:1) [24]. Even in a gas with 
fixed chemical composition, »vressure and translational tem- 
perature, and hence, fixed Ap/p, An/n can vary according to 
the ionic and electronic excitation of the medium caused by 
pumping. 

Consider the Ladenburg formulation of the Kramers-Kronig 


dispersion relation [25] 





ak N, i x NS. 
<= — a = Pe nares 

Er i aD Dele meee a et) 7 Spee! 

™m™ Cc i ee age 

i] 2 k kj 
(3) 

e electron charge 
‘a speed of light 
m mass of electron 


total number 
See wavelength for transition between k and j 
A wavelength of interest at which n is derived 
ae MomOcc@llatOr Stmeneth for transition from k to j 
is degeneracy of ees energy level 


N. population in een energy level 


There are 2 distinct species indicating specific electronic 
or ionic excitation present in the cavity. Index } sums 
over lower levels and k over upper levels of possible con- 


figurations. Identify the mole fraction as 
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the fraction of the total molecular population represented 


by species 2. Finally, define a quantity 


re 


oe ia Nh Si 
‘) 3 Reef [1 + (7) 
ai 1 ; A = dD. nes 
e j Ks kj 
In terms of the quantities defined above 
n-l=N )) X,T, (8) 
i 


Seeeene Index of refraction in the cavity is the sum of the 
contributions of all the species present in the cavity. 


Since Ny - NX, 


oTy + ... (9) 


The partial density is given by Pies: Nom) where mM, is the 


molecular weight of ae species. In terms of partial 


density, equation (9) becomes 


n-l = 0,7, /m, + Py t,/m, ae CG) 


Introduce the mass fraction for the gas mixture Yo = p,/0, 


where p is the density of the gas mixture. The result is 


n-l = o[¥,T,/m, + YoT5/My | (Giei) 


Le 





Equation (11) shows that n-l is proportional to 0; however, 
equation (11) further demonstrates that n-1l depends on mass 
fraction and index of refraction of each species. Even with 
uniform flow, the spatial phase relation across the cavity 
will not be uniform. Accurate modeling of the phase effects 
requires knowledge of the species populations and the con- 
tribution to the index of refraction of the species in the 
cavity. This thesis will explore a model for determining 
the species populations ina laser cavity and apply a model 
originally developed at the Naval Research Laboratory (NRL) 
to calculate the populations of the different configurations 
of excited Neon present in a XeF laser. 

Borst, in section IL, a review of electronic transition 
gas laser kinetics is discussed. Second, in section III, a 
simple model for an excimer (KrF) laser is developed. Third, 
in section IV, the computer program developed and written at 
NRL by Johnson, Palumbo, and Hunter [26] is discussed. 
merch, in section V, a formulation of input reactions to 
the NRL program to calculate the populations of excited 
meen Configurations is developed. Fifth, in section VI, 
@eeculation of the index of refraction is performed. Sixth, 
in section VII, the calculation of phase shift is performed. 


Mmieety, in section VIII, the results are discussed. 
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Electronic transition gas lasers (i.e. exciplex or 
excimer) generally involve complex, branching chains of 
energy flow from initial pump excitation of background gas 
to eventual extraction of laser radiation. Figure 2 illus- 
trates major energy flow pathways in e-beam pumped 
Ne/Xe/NF. mixtures. These energy chains include processes 
such as: energy deposition by electron ionization and 
excitation; energy transfer by various collisional 
mechanisms, including neutralization, displacement and 
“harpooning", leading to formation of the upper laser 
level; andcollisional mechanisms which intercept excitation 
energy before it can reach the upper laser level, or which 
quench that level before stimulated emission occurs. Also, 
one must consider extraction of stimulated emission includ- 
ing absorptions by various processes in the gas and reso- 
nant reabsorption by the lower laser level (if there is 
one), laser cavity oscillation, and coupling through output 
mirrors. A time-dependent model of the kinetics of such a 
system must define and update species densities for a large 
number--typically between 20 and 69--gas components, includ- 
ing ground states, several excited states and ions of atoms 
and molecules, plus electrons and photons. 

The time evolution of these components is followed by 


me solution of a set of coupled first order ordinary 
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Figure 2. Major energy flow pathways in e-beam pumped 


Ne/Xe/NF3 mixtures (Reproduced from Huestis [27]). 
Note: Energy input from the electrons comes in along 
the pump arrows. The species in each box, plus 
the comvonent along the reaction arrow, yield 
the result in the next box. The wavelengths of 
the emitting species are indicated. They may be 
Cie Cemeositenentmne by Mr or Xe (not shown) . 
NeF* is presumed to predissociate. 





differential equations; one equation is required for each 
gas component. To the set of equations defining species 
densities are added further equations describing the evolu- 
tion of the photon field, the gas temperature, and (in some 
models) coupled electrical circuit parameters. The species 
equations can be generally described as equating the rate 
of change of species density to a sum of contributing for- 


mation terms less a sum of depletion terms as 


cae 2 Pui” = Date (12) 
ic 


Here N. is a particular species density, the Pai are the 
formation terms, and the Dey the depletion terms. For 
example, consider a simple chemical reaction of the form 
A+B-> C+D whose forward rate constant is k: if the 
Species densities are expressed in om, the units of k are 
om?/sec. For this reaction, the rate equations defining 


species C and D would contain the formation term 
F = k{A][B] ; CES) 


at the same time, the equations defining species A and B 


would contain the depletion term 
D= k{A][B] . (14) 


Recall that square brackets for a chemical symbols means 


number density of the chemical with dimensions L/om>. 


ZA 





Consider equation (12) again. Species iis formed in 
j different reactions. Likewise, species i is depleted in 


k different reactions. 


AGE 





ipo eee INET IGS MODEL FOR KrF 


In a thorough model of an excimer laser a large number 
of kinetic reactions are necessary; see, for instance, the 
Momereactions listed in the output of the NRL program for 
XeF (Appendix A). The reactions are used to obtain the set 
of coupled first order differential equations. The solution 
to this set of equations is very involved and tedious 
because of the difficulty in handling the tangle of reactions. 
Therefore, only a few of the kinetic reactions in a KrF 
laser are considered. The object of the simple computer 
model was to obtain a qualitative insight to excimer laser 
finetics. 

The model uses a He/Kr/F, etemre megs! / 1. 77/0.13) with 
He as diluent. Figure 3 illustrates a simplified e-beam 


pumped KrF laser geometry. 


ELeCTRON BEAM | 


ae 





+ — 


Point of Stimulated 
Enission 


LASER CAVITY 


be cars LENGTH (L) ————> 


Figure 3. Example of e-beam pumped KrF laser geometry. 
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The first step in the kinetics chain leading to extrac- 
tion of laser radiation is the deposition of energy in the 
diluent gas by either photons or electrons. In this model 
electrons are used. The electrons'‘ionize the Kr and 


dissociate the F, 


k 
e' + Kr +t Kee + e+ e' Gir!) 


where e' is a primary electron, e is a secondary electron 
resulting from ionization of Kr, and ky Womene. fLoLward 
rate constant. The secondary electron has sufficient 
energy to dissociate and electronegatively attach to 
eeomic F. 

s 


ese E+ F- (16) 


iteeaacdition dissociation of F., by the primary electrons 
Seecurs. 


Ka 


Having produced the ion pair states, recombination of 
positive rare gas ions and readily formed halide ions yields 


the excited state (KrF*). 


Ky 
v . 
Kr +F +™M => KrF* + M (18) 


The KrF* can produce a photon by stimulated emission. 


k 
KrF* + hy => KrF + 2hv (19) 
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Pubs@, Lecombination of F can occur 


Oye 


Je ieee a eg 1 eae a (20) 


Fluorescence decreases KrF* population by 


E 
KrF*¥ + KrF + hv (21) 


Where t is the radiative lifetime. 

Reactions (15) to (21) are the only reactions considered 
in this model. None of the quenching processes of KrF* 
have been considered, such as 


k 
KrFe + 2Kr > KraF* + Kr Ce) 


Z 
nor have many other reactions [2] and [28]. 

The reactions considered in this model are used to find 
the formation terms equation (13) and depletion terms 
equation (14). For example, examining reaction (15) the 
formation term for the concentration of e and Kr’ is 


@erermined. 


Formation e k, le’ ](Kr] (23) 


Formation Kr’ k,le° ] [Kr] 2) 


where [ ] denote concentration in specie number/cm>. 
The depletion term for concentration of e from reaction (16) 


is 


8, 





Depletion e = ko lel [Fo] >) 
and for Kre from reaction (18) 
+ +, 0- 
Depletion Kr = k, [Kr Woe 80) (26) 


Applying equation (12), that is, equating the rate of change 
See concentration to the sum of contributing formation 
terms equation (23) less the sum of depletion terms equation 


(25) yields 





ciel = kyle] {¥%r] - ky lel [Fy] (27) 
Pummeariy, combining equations (24) and (26) 

d{Kr’] ee 

G1Kr ] = fet] [Kr] - k,[ke*] [F] [0] (28) 


Wemeinuing the process 








SEL = kylel[F,] + kgle'I[F,] - kg(FIIM] (29) 
ere ale ith.) - &, (kc CF 1M) 
dt 2 2 3 - 4 
(30) 
ctRePel = k,fKeT](F JIM] - ko hv] [KeF*] - [KrP*]/+ 
(31) 
SLKrF = ke [KrF*][hv] + [KrF*]/t eZ) 


are obtained. 
In reaction (15) the e' is the e-beam input and there- 


meme 1S a function of time. The model assumes the 
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concentration of e' is initially 1.9E13 electrons/em> then 
becomes zero after about 177 nanoseconds. The secondary 
electron mean energy is 10 eV. 

The rate constants (k) for the reactions were taken 
meom seclentific articles [2, 29, and 30] or derived as in 
the case of Ke. The values of the rate constants are listed 
in Table IIL. The radiative lifetime, t, is equal to 15 


nanoseconds from reference [2]. 


imope Li, Rate Constants 


Rate Constant Value Reference 
k, 5.0E-09 em?/sec 29 
k. 2.0E-09 em?/sec 30 
k. 5 .0E-10 i BEE 30 
ky 1.0E-25 cm°/sec 2 
ke 1.5E-05 ome fae esbaabe pee 
ke 1.0E-31 em?/sec 2 | 


Equations (27) to (32) can be solved to calculate the 
concentrations as a function of time. Multiplying these 
equations by dt yields a form for numerical solution on the 


Me2o50. Equations (27) to (32) become 


dfe] = (k,[e']([Xr] - kyle][Fo])dt (33) 


ae = 


aKr | ies | [FF iieedte (34) 


(kale'j(Kr]) - k 


1 4 


Zo 





d(F] = (kjlel(F,] + k,le'][Fy] - kg(F][M])dt (35) 
dIF°] = (kylel[F,] + kale'}(F5] - k,[Ke ](F-][M]) de 
eG) 
d[KrF*] = (k, [Ke ][F ](M] - ke[hv][KrF*] - [KxF*]/t)dt 
(37) 
d{KrF] = (ke [KrF*] [hv] + [KrF*]/t)dt (38) 
respectively. 


The initial conditions for the set of equations are: 


[e] = 0 (39) 
[F] = 0 (40) 
[12 cea © (41) 
[kr] = 0 (42) 
fr} = "0 (43) 
[KrF*] = 0 (44) 


Several constants are specified as follows: 


[hv] = 4.0E11 (calculated in Appendix C) (45) 

[F.] = 4.0E16 (46) 
Assumed values for 

[Kr] = 5.4E17 } the He/Kr/F9 mixture (47) 
Cet 7/ 0-13) at 

[He] = 3.0E19]} STP. (48) 


fimitial T = 0 


final T 250 nanoseconds 


Delta T =meminal © > Inrtial 1) /256 

A computer program was written using equations (33) to 
(38) and the initial conditions and appropriate constants 
given by equations (39) to (48). Figures 4 and 5, which are 


plots obtained from the HP9$30 program, are the concentrations 
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in specie number/cm> versus time from 0 to 259 nanoseconds. 
The plots in Figures 4 and 5(b) indicate that some species 
densities rise extremely steeply in the first tenth of the 
simulation time; later in the pulse (after about 50 nano- 
seconds) species densities tend to flatten out. See Figure 4. 
Once the e-beam terminates, the decay time is approximately 
50 nanoseconds. In Figure 5(a) and (b) other species 
densities increase gradually and flatten out. The results 
of this model are in qualitative agreement with the NRL 
Simulation contained in reference [26]. The NRL simulation 
limited the change in concentration to 5 percent per time- 
Seepeco Obtain sufficient accuracy ({26]. The question of 
Smesteps for the simple model was avoided; sufficiently 
small steps were used throughout the program. 

Table III lists the variables, symbols, and equations 
used in HP9830 computer program. 

Table IV lists the KrF Computer Program and Table V 


Mees the output. 
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Kegu molecules/cm>, as a function of time. 
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Table IIil. 


HP9830 


Kl 
K2 
K3 
K4 
K5 
K6 
38 


Variables, Symbols, and Equations Used in 


HP9830 Computer Program for KrF Laser. 


Variables and Symbols Used in the Program 


Text 


e 


1) 


Meaning or Value (for constants) 


primary electron 

secondary electron 

Fluorine atom 

Fluorine molecule 

Fluorine atom with extra electron 
Krypton 

Krypton minus an electron 

sElnaligval jeyevel ty 

Krypton Fluoride 

Krypton Fluoride excited 

photon 

Rate constant 

Rate constant 
Rate constant 
Rate constant 
Rate constant 
Rate constant 

Radiative lifetime 

I=1 for initial conditions 


Index; 


Time increment 


lige tes 


joules 

om?/sec 
aE eee 
cm/sec 
ant figac 
em?/sec 
em’ /sec 


Sha 


sec 





Table III (Continued). Variables, Symbols, and Equations 
Used in HP9330 Computer Program for XrF Laser. 


Equations Used in the Program 


Line Number Equation Number 
860 58 
870 34 
880 55 
890 36 
900 45 
910 Ws 
B20 38 
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IV. XeF LASER COMPUTER PROGRAM FROM TRL 


In order to calculate the excited states of neon in a XeF 
laser a complex kinetics program is required. A kinetics 
scheme for the XeF laser developed at the Naval Research Lab 
Peamused. it was written by Dr. Louis J. Palumbo [26 and 31] 
who assisted greatly in modifying it for use on the IBM 360/67. 

The computer model consists of two parts: first, an e- 
beam deposition section, in which the kinetic energy of near 
relativistic electrons is used to ionize and excite the laser 
gas, and, second, a chemical kinetics section, in which the 
energy deposited in the gas cascades down by various processes 
to the upper laser level and is extracted as laser radiation. 

One input to the NRL computer model is a list of reactions, 
Seen the format of equations (16) to (21). Also an input 
is a list of reaction rate coefficients. The computer program 
sorts the reactions and forms formation terms, F, and 
depletion terms, D. Formation and depletion terms were 
discussed earlier in connection with equations (13) and (14). 
Continuing, the computer program selects the appropriate F 
and D terms to form the rate equation of equation (12). Next, 
the computer program multiplies the rate equation by a variable 
timestep. The same approach was used in the simple computer 


preeetam; see equations (33) to (38). 
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An advantage of the NRL program is the ability to change 
the number of reactions. The computer program does not change; 
only the input changes. The differential equations are 
generated internal to the program. This method makes it much 
easier to alter the set of equations and greatly reduces errors 
attendant to such alterations which arise with ease in the 
complex tangle of reactions. The detailed explanation of the 
computer model is contained in references [26] and [31]. 
Appendix A gives a list of reactions used to generate the 


population distribution within neon energy levels. 


~~ 
fee 





V. CALCULATION OF NEON EXCITED 


STATE POPULATIONS 


The excited states of neon are represented in Figure 6. 
The electronic configurations of interest to this study are 
mtemos, 3p, 4s, 3d, 4p, 5s, 4d, 4f£, 5p, 5d, and 6s. 

In order to calculate the populations of these particular 


configurations, one needs to know the reactions involved and 


the rate constants. For example, 
al 
Ne + e-beam > Ne(3s) 
o 


Ne + e-beam + Ne(3p) 


k 


irae 
Ne + e-beam + Net5s) 


These reactions could then be input to the NRL XeF pro- 
gram to calculate the populations of the excited configura- 
tions of neon. After, an extensive literature search (Optics 
and Spectroscopy, The Journal of Chemical Physics, Applied 
Physics Letters, Physical Review A, Journal of the Optical 
Society of America, Chemical Abstracts, IEEE Journal of 
Quantum Electronics, etc.) rate constants for the individual 
Seemted configurations i.e. 3s, 3p, ... were not found. 

Therefore, it was necessary to make a number of assump- 


tions. The procedure followed was to extract all reactions 
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Figure 6. Energy-Level diagram of Neon (Reproduced from 
Pressley [32]. 
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from the NRL program involving the neon excited state (Ne*). 
Then, these reactions were subdivided into two groups, the 
Ne(3s) and all higher configurations, which are denoted by 
NE*1 and NE*Z respectively. The population distribution 
within NE*Z was calculated assuming a Boltzmann distribution. 
The list of reactions which were used as modified input 


to the NRL program are listed in Table VI in the following 


format: 
Rate Constant 
Ne + e' + Ne*¥ + e' I.20E-18 ~~ (49) 
Ne + e' + Ne*l + e' 0.925E-18 (99) 
Ne + e' + Ne*Z + e' O7275E-Ls (51) 


To compare equations (49) to (51) with the first three 
reactions of Table VI note that e' is ZBEM in FORTRAN nota- 
tion. The other symbols are obvious. The first reaction 
yielding Ne* is divided into the Ne*l and Ne*Z reactions. 

A computer run with the current density equal to 10 A/cm? 
and a computer run with the current density equal to ll 
A/ cin? calculated the populations for Ne*l (3s) and Ne*Z 
(all levels greater than 3s). 

The rate constants for the reactions were determined by 
taking a fraction of the total rate constant. For example, 
in table VI, for the second set of reactions, the rate 
constant for NE*l is 0.8 and NE*Z is 0.2 of the rate con- 
Stant for NE*. These fractions of the rate constant for 
NE* are noted in parentheses following the reactions. The 


choices of the fractional values were arbitrary. 
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“ 
TABLES VI. Input Reactions to the NRL Program. 
Reaction Rate Constant 

NE + EBEM + NE* 52 SI iO) = 1.8 
NE + EBEM > NE*1 EM 0.925E-18 
NE ae eyeyl re ees Z, eo! 0.2752-18 
NE* + E- ene + E- 2.00 E-10 
Meet 2+ CL sa |N0. + E- (.3) 1.60 E-10 
ies, +t i- See NE, + E- Cez,) 0.40 E-10 
NE2+ + E- > NE* ae Ue Sey a0 2105) 
NE2+ + E- eo: NES Cao) 3.15 E-08 
NE2+ + E- > NE*Z a2 Jes CT) 0.35 E-08 
NE3+ + E- > NE* a YE az Nid 3.50 E-08 
NE3+ + E- a2. A(ideaal ae + NE G5) 3. 3252-06 
NE3+ + E- aN a Z sa a= NIE C05) "0. 175-03 
XENE+ + E- > NE* a DID 6.00 E-08 
XENE+ + E- > WE] ae e C0) S00 20'S 
XENE+ + E- + NE*Z ina G0 (a4) 2.40 E-08 
NEX* one, + NE > NE2* a= S05; 520P eso 
eer] + NE + WE > NE2* stun) op 4.0 E-34 
fer «2t NE + NE > NE2* + NE Cre) PO ee 
NE* Sek = ao + E- se Ne 7.40 E-1ll 
Mex) + XE =~ XE+ + E- ee, (.6) Cpe ce y= 
NE*XZ + XE == ).J2ae + E- + NE Cua) 2.96 E-ll 
NE* ee E + NEXE+ + E- DO: = 
NE*] + 2% > NEXE - B= (.6) JL6 SVAee oye 
NE*Z + XE + NEXE+ + E- C4) 0.88 E-1l 





teases Vi (Continued). 





ie ttenedacclons to the [RL Program, 


Reaction Rate Constant 
NE* + F2 +> NEF* + F 4.10 E-LO 
NE*lL + F2 +> NEF* + F (C23)) Bao = 0 
Me<. + F2 +> NEF* + F Ge2) 0.62 E-L10 
NE* + NF3 > NEF* + NF2 1.05 E-LO 
NE*L + NF3 > NEF* + NF2 Casp 0.84 E-10 
Nes. + NEF3 > NEF* + NF2 (eZ) Oree2 ik B10 
NE* eee 2 > NEF* + NF 1.0 E-LO 
NE*] + NF2 > NEF* + NF CaS) 0.8 E-10 
NE*Z + NF2 > NEF* + NF a2) 0.2 E-10 
NE* + NF > NEF* + N 1.0 E-10O 
Mee) + NF > NEF* + N (.8) 0.8 E-LO 
Ne~Z + NE > NEF* + ? (a) O22) E10 
NE2Z* + HN352L > NE*¥ + NE + ABSL L.0 E-l18P 
NE2* + HN352L > NE* + NE + ABSL erg ) 0.9 E-18P 
Mee~ + HN352L -> NE* + NE + ABSL Col) O.l1 E-18P 


* 
Use of capital 
eeimtine. 


letters more closely represents computer 
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However, the first reactions listed in Table VI do not 
have fractional values in parentheses, as the reaction rate 
constants were determined using reference [9]. The plot of 
"Excitation of Neon by Electrons" from page 467 of refer- 
ence [9] was extrapolated to 10 keV, and the value for the 
2P1 state was assumed to represent all excited configura- 
tions above the 3s (NE*1). The difference between the rate 
constant for NE* and the 2P4 value yielded the rate constant 


for NE*1, 


i 


—_—- ee wae ee ree eee 





Vien ecurattON OF THE [INDEX OF REFRACTION 


The values (n-1)/N where n is the index of refraction 
and N is the species population, have been calculated Refs. 
me (34) and [35]. Figure 7 lists the values of (n-1)/N. 
References [34] and [35] contain the details of calculating 
these values. 


The index of refraction (n-l) is equal to 





_ ,(n-1) n-l: n-l 59 
ground ¥* * 
ae Ne] Ne*Z 


Using the populations calculated in the NRL model for N, 
N(3s)., and N(Z) equation (52) can be solved. However, the 
populations of individual electronic configurations must be 
determined. 

One feature of excimer lasers, in contrast to the 
vibrational transition lasers of the infrared region, is 
the extent of randomness of populations. In the cO, gas 
dynamic laser, mode temperatures can be assigned to the 
various vibrational modes. Due to the short span of excita- 
tion and power extraction in an excimer laser, the various 
modes do not have an opportunity to equilibrate. This fact 
is the reason extensive computer programs are necessary for 
kinetics. An assumption of a Boltzmann distribution is 


unfounded. However, in view of the sparsity of reaction 
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rate coefficients for excitation to specific configurations 
a Boltzmann distribution is assumed. Assuming a Boltzmann 
relation for the upper levels the populations of the 3p, 


Pomoad, ... 6S can be calculated. 





ee — Nt3p) + N@4s) + N(3d) + ... + N(6s) Ge) 
~ (Ey57E 3p) 
Bhs — 
N(Z) = N(3p) [1 + — exp 
8 3p 
‘ ~ E3q7F 35) 
eeclexp 9 eT C4) 
3p 


where g denotes the degeneracy of the configuration which is 
equal to the sum of the 2J + 1 terms. E is the energy sub- 
memeoeed by the configuration. 


pelving for N(Zp) 


WCE, 
N(3p).= === (35) 
4s ~3p 
Ll + ~ 4s exp ki a 


All the quantities on the right side of equation (55) 
mem <nOwn Or can be calculated. Table VII lists the values 


of the degeneracies and energies for the configurations. 


Die 











TABLE VII. Degeneracy and energy values for the 


configurations. 
Electronic 
Configuration 
3s eZ 
3p 36 
4s WZ 
3d 3 
4p 36 
5s eZ 
4d Ve 
4f£ 84 
3p 36 
5d 18 
6s eZ 


Remes +l) 


134674. 
EOI 1s. 
I oysie) 0) ¢ 
HOLS 7G. 
bo3lso. 
ONG a 


166500 
Ooo 


166 Joie 


E (om + 


168000. 


168750 


The temperature T is calculated as follows using the 


Boltzmann relation 


= (080 SJ, 
SS 20 
N(3s) _ 53s fal kT 
N (2p) Bop 
Semying for T 
a, St a 
N(3s) 82p, _ ~S#357E 2p 
2a NCP) 34_ kT 


(56) 





= (3g) ape, 
T = __ Se ee ee (57) 


N(3s) 89 
k &n Np) g,.! 


Siiostrtuting the value of T into equation (55) yields a 


value for N(3p). To calculate N(4s) one uses 


g Es 30" 
N(4s) = N(3p) —% exp kT (58) 


Similarly the populations of the other configurations 
Can ve calculated i.e. lI(3d), N(4n), ete. Once the povrula- 
tions o£ the configurations have been found it is a simple 


calculation to obtain n-l applying equation (52) as 


n-1 n-L n-L 
Beep)? tle(4s5) * °°2 He(Ss) 


are given along the bottom of Figure 7. 

A computer program to calculate the index of refraction 
(n-1) using equations (52), (55), (57), and (58) was 
written. 

Pmeorsiee) / the 35 + 94p transition contains an R*. R* 
indicates that resonance may occur and can have a larze 
effect on the index of refraction. The spikes in Figure 8 
are located at the resonant wavelengths. Resonance is 


mecsewssed in reference [36]. 
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The plots in Figure 9 are of the calculated indices of 
refraction (n-l1) versus time (nanoseconds). Figure 9{a) 
is tor a laser wavelength away from 3s + 4p (no resonance) 
for both low and high power. The current density for low 


2 and the hizh power is ll A/cm’. For a 


power is 10 A/cm 
nonresonant wavelength, the change in n-l is in the fourth 
decimal place. 

The resonance effect was taken into account in a second 
computer run; see Figure 9(b). For a resonant laser wave- 
length, the contribution to the index of refraction due to 
Ne(3s) was calculated assuming (n-1)/Ne(3s) equal to 
-5.0208; the corresponding value for Ne(4p) was (n-1)/Ne(4p) 
epee O | .66/E9G. Figure 9(5) is for a laser wavelength 
near the 3s + 4p resonance for low and high power. The 
resonance has a marked effect on the index of refraction 
in comparison to Figure 9(a) as evidenced by the ordinates 
of the plots. Also, Figure 9 reveals that increasing the 
power decreases the index of refraction. 

The right hand ordinate scale for Figure 9(a) is the 
iat 10 Gary m1). ; where Goa 1)) is the value at zero time. 
The scale emphasizes the fact that the relative changes are 
in the fourth decimal. 

Table VIII lists the variables, symbols, and equations 
used in the H?P9830 program to calculate n-l. Table IX lists 
the 4P93830 program used to calculate n-l and Tables &% and XI 


list sample output for low and high power respectively. 
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(a) Index of refraction (n-l) at laser wavelength away from Ae + 4p 
resonance for high and low power. 
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(6) Index of refraction (n-l) at laser wavelength near 3s + 4p resonance 


for high and low power. 


Figure 9. Calculated index of refraction (n-l) as a function 


of time. 
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TABL™ VIII. Variables, Symbols, and Equations Used in 
PE os Ueeempucer =rorram for Calculation of 
n-l. 


Variables and Symbols Used in the Program 
HP93 30 esc Meaning or Value (for constants) ines 


G g degeneracy (2 2J+1) = 
E E energy em” > 
Gb a speed of light cm/sec 
K2 k Boltzmann constant eV /°K 
M n-1/Ne() COMEalDOuUELOM OL excited state - 
K BaD geo Sa sis Index of the levels 1 to 12 - 
Ke) aE Temperature degrees 
Xelvin 
J - Time increment index - 
Hqguations Used in the Program 
Line Number sguation Number 

320 Dy 

730 56 

O'30 54 

1050 55 

1130 53 

1359 D2 
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TABLES X. Output of the HP9830 for n-1l at Low Power. 
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TABLE NI. Output of the HP9830 for n-1l at High Power. 
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Vibes Cec uULATION OF PHASE SHIFT 


Beam quality as previously mentioned depends on favor- 
able illumination of the exit aperture; favorable illumina- 
tion implies near constant phase. Therefore, calculation 
of the phase shift is important. To calculate the phase 


shift consider a laser cavity of length L as shown in Figure 


OF 


CAVITY 
LASER OUTPUT 
 .——— 
Figure 10. Geometry of Cavity 


Other symbols essential for the calculation are 


N - number of waves AX - wavelength 
1 - denotes low power n - index of refraction 
2 - denotes high power vy - frequency of the laser 


The number of waves in length L is given by 


N, = L/A, (59) 
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iemestaeion (22) i = 4,2. Wavelength is related to index 


Se rerraction by 


Ay = Ways ae (60) 


Consider two regions in the cavity which are pumped at 
different power inputs. The difference in the number of 


waves in the two rezions is 


No - M, = LG - 
1 


) (61) 
2 a No 


Substituting equation (60) into equation (61) yields 


ry) 


Mo - Ny = L(— - — = = (,-n)) (62) 





Multiplying the right hand side of equation (62) above and 


below by n, produces 


(63) 


Multiplying equation (93) by 1/L and substituting equation 


(60) yields 


No = Ny a No-Ny 
L X 





) (64) 


ee 
2) 
tr 


The term on the left hand side is the phase shift per unit 
length of cavity. The phase shift at 1000 nanoseconds for 


a nmonresonance wavelength is now calculated. 


} = 353 nm = 353.0E-7 cm (Figure 8) 





Beom Dahbles = and XL the values of n are 


ny = 1+ 5.420392-05 


ny et 5. 42049 e=05 


Substituting into equation (64) 


Noh __1___(1.9000542039-1..0900542049) 
L SoU = 7 MOO W05 42049 


= 29 6207923 Bae 


The phase shift at 1900 nanoseconds for resonance is calcu- 
lated below. 


ie oo) ty Ob 7 


mae values of n are 


no Pee E7306) 


ny Perec 025-06) ) 


Substituting into equation (64) 


opel a i Po aemeiaoy LO le. JOC00560Z 


1 
ames  35n.cn-7 *—  1.000006602 


y= -3.959E-3 cn 


The Strehl ratio, which is discussed by Born and Wolfe 
[21], provides a criterion for beam quality in terms of 


phase distortion. The ratio is 


I= 1, exp (-(216/a)*] (65) 
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where I is beam intensity in the far field due to a phase 


distortion of 6/A. Based on equation (65), the criterion 











et 
o/s 0.1 (66) 
is frequently invoked for satisfactory beam quality. Com- 
paring equations (65) and (61) one notes that 
M5 - Ny = 6/d Co 
The waximum lenzth of the laser cavity, L, for good 
beain quality is | 
(N5-N, ) 
N-N NA au 
Ve - 
(4) 4S 
Substitution of numerical values yields the following 


results: 
- maximum L for a laser operating at a nonresonant wave- 
length is 35 meters. 
- maximum L for a laser operating at a resonant wave- 
Pength is 0.1116 meters. 
For a nonresonant laser wavelength, a laser cavity with 
length, L, less than 35 meters will have 6/A less than 0.1. 
However, if the laser wavelength is in resonance with neon, 
a length of laser cavitv of only 0.1116 meters will cause 
OrrL: 


6/A = Recall the power into the laser cavity was 


Varied by 10 percent. 
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Vil te BoUlG S 


Yor a power variation of 10 percent the phase shift at 
nonresonance is small (No - Ny = -9.0028 for a 199) cm 
cavity) and is dominated by the neon ground state whose 
population is approximately one million times larger than 
the other species. However, in the resonance case, if the 
assumed values of R* are correct, the index of refraction 
can be greatly magniried, resulting in a phase shift of 
morgage 3 em) as compared to -2.8327E-5 en + for nonreso- 
nance. 

It should be noted that the populations of the neon 
excited states were calculated using an arbitrary choice to 
determine the rate constants. However, the population of 
Ne is approximately six orders of magnitude larger than 
the neon excited states and thus dominates. As long as the 
laser is operated away trom the 3s > 4p resonant wavelengths, 
the phase shift will be negligibly small resulting in 


Satisfactory beam quality. 
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CALCULATION OF RATE 


APPENDIA B 


CONSTANT Ks 


Consider a thin layer of active optical medium, with 


=a 


gain a, cm 


change in intensity I, 


, in the laser cavity of length 2, 


cm. a2he 


watts/cm¢, within the thin layer is 


(09) 


The number of photons per unit time entering the layer of 


active optical medium is 


° ieee" 


yg 


re 


A corresponding term for 


Mimic: time is 


out hv 


The rate of production 


the difference between 
{rate of pro- : 
me tion = Se 
lof photons | 


Equation (72) neglects 


. — tout 


(79) 


the number of photons leaving per 


Ce 


of photons within the thin layer is 


are and Ne which is 

ao 7 me a (I = ale ) (72) 
ng ge We in : 

spontaneous emission. Combining 


equations (69) and (72) yields 


76 





production rate = a (73) 


From equation (19) the rate of production of photons by 


Sermitated emission is 


production rate = k,[hv][KrF*]Ag GC 


where the symbol [hv] means number of photons per unit volume. 


The gain coefficient also is equal to 
a = o[KrF*] Ge) 


where o is the optical cross section for stimulated emission 
at the frequency v. The number density of photons [hv] is 


related to intensity by 


eon) = = (76) 


mere c 1S the velocity of the photons; c, obviously, is the 
Speed of an electromagnetic wave. 


Combining equations (73) to (/6) gives 


3 
meee dimensions of k. are cm /sec; the dimension of o is cm, 


D 
and cis ei eee From reference [2] the value of o is 
5. 055-16 om, eaceceis 3.9H10 em/sec. Substituting these 


values into equation (77) yields 


al 





ke = 1.5E-05 erin cee 


This is the value for Kk. used in the computer program. 


7€ 





NNSA eas BBR Cs @ 
CALCULATION OF PHOTON CONCENTRATION 


The concentration of photons is given by equation (/6) 


as 
[hv] = L/chv 
Assuming, I = 10 any ents (saturation value) 
c = 3.0H10 cm/sec 
h = 6.6256E-34 J-sec 
Me = 243.0E-9 m 


the concentration of photons can be calculated. 


hv] LORoES 


(3.0E19) (6.6256E-34) (spa*pee) 


=e soniivicn- 
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